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1.  INTRODUCTION 


One  class  of  problems  of  current  military  interest  concerns  the  * 

response  of  antennas  mounted  on  the  exterior  surface  of  metallic  ^ 

equipment  shelters  to  the  tactical  electromagnetic  pulse  (EMP) 

j_. 

environipent  created  by  a near-surface  nuclear  detonation.  With  the  hope 

ij 

of  ultimately  obtaining  models  suitable  for  engineering  purposes, 
walytic  formulations  of  such  problems  are  posed  in  ideal  geometries 
with  simplified  problem  physics. 


One  such  formulation  consists  of  a perfectly  conducting  sphere  of 
radius  a,  in  unbounded  space  filled  with  a medium  of  constemt 
conductivity  o,  in  the  presence  of  Compton  electron  flux  of  velocity  v 
euid  current  density  J^.  The  source- current  density  is  curl-free  and  the 
total  electromagnetic  current  density  (in  the  absence  of  the  sphere) 
vanishes.  The  sphere  is  treated  as  ideally  gamma-thin;  that  is,  for 
each  electron  entering  the  sphere  another  electron  leaves  the  sphere  in 
the  direction  of  the  photon  flux. 


A formal,  frequency-domain  series  solution  to  the  latter  formulation 
has  been  obtained  by  Bombardt.^  In  this  report,  we  obtain  an 
approximate  time-domain  solution  for  the  electromagnetic  fields  on  or 
about  the  sphere  by  contour  integration  of  Inverse  Laplace  transforms. 
The  approximate  solutions  are  used  to  perform  some  tentative  estimates 
of  the  coupling  of  the  radial  electric  field  near  the  sphere  to  zm 
emtenna. 


2.  SPHERE  MODEL  AND  APPROXIMATE  SOLUTION 

The  problem  geometry^  for  the  sphere  of  radius  a is  shown  in 
figure  1.  We  use  mks  units  with  free~space  permittivity 
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(fo./ifl.O) 


Figure  1.  Sphere  exposed  to  Con^ton  current  density — A/m^  in  presence 
of  constant  conductivity  o . 


S,  Bombardtf  Jr.,  Electromagnetic  Response  of  a Sphere  Over  a 
Ground  Plane  in  the  Presence  of  Source  and  Conduction  Currents  in  the 
Air,  Barry  Diamond  Laboratories  TR-1716  (October  1975). 
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e - 10“®/36it  F/bi  and  free-space  penneedilllty  ii  ■ 4Tr  x lo”^  H/m.  The 
o o 

Cc»^ton  current  density  is,  with  explicit  sign  for  the  electron  charge. 


-^c('^  - f)  " 6 - f)  ^ 


- f(a))e^®*,  6 


u 

V 


(1) 


u(t)  - 0 (t  < 0);  = 1 (t  > 0)  . 


In  equation  (1) , 
t =>  time, 

z > horizontal  coordinate  (see  fig.  1) , 
u ^ Heaviside  unit  step, 
u = radian  frequency,  and 


-f((o)  - /*j  (t)e^“^  dt 
— »*'c 

-j^(t)  = /*f(u>)e’^“^  du) 

c 2ir  -«» 


(2) 


The  electric  field  in  the  medium  in  the  absence  of  the  sphere  is 

iBz 


= 

z 


•f  (u))e 


a - iue 


(3) 


The  formal  series  solution  of  the  boundary  value  problem^  yields  the  sum 
of  incident  euid  scattered  electric  fields: 


total 


^i  . s 
= E + E 


(4) 


where 

H {r,0  ,u) 
9 


-f(u)  y 

V 


n-1 


(2n  + 


— n»l 


(5) 


^J.  W.  Bonibardt,  Jr.,  Electromagnetic  Response  of  a Sphere  Over  a 
Ground  Plane  in  the  Presence  of  Source  and  Conduction  Currents  in  the 
Air,  Harry  Diamond  Laboratories  TR-1716  (October  1975). 


t 


I 

) 

i 

?■ 


I 

I 


^ ! 


. i 


fi- 

i 


, I 


I 


i 


i i 

I 


Lr 

[Lt-  ‘ 


E“(r,0,(i>) 
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f(a)) 


- S i""^n{n  + 1) (2n  + 


(6) 


IrK"  ‘“>] 


u«ka 


Eg (r,0,u) 


f(<«>) 


.n-1 


(2n  + l)j 


n(^)- 


du 

uh^^' (u) 
n 

u=kr 

du 

uh^^^  (u) 
n 

u«ka 

(7) 


3p 
n 

■5^ 


vrtiere  j is  the  spherical  Bessel  function,  P the  Legendre  polynomial 

(1)  " 
and  h'  ' 


jj  jjj  + the  Hankel  function  of  the  first  kind.  All  special 

functions  are  defined  with  the  conventions  of  Abramowitz  and 
2 


Stegun.  For  convenience,  we  convert  the  Fourier  forms  in  equations  (5) 
to  (7)  to  Laplace  form  emd  perform  (term  by  term)  an  evaluation  of  the 
inversion  integral 


d+ioB 


g(t)  = 


5r  / 

d-ioo 


(8) 


where  9(2^)  represents  the  desired  field  from  equations  (5)  to  (7)  and 
Zq  is  the  Laplace  transform  variable,  s,  regarded  as  a complex  varledsle. 
A Bromwich  contour^  is  employed  (see  sect.  5). 


Abramowitz  and  I.  A.  Stegun,  eds..  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards,  AMS  55  (1965). 


F.  B.  Hildebrand,  Advanced  Calculus  for  Engineers , Prentice  Hall 


(1949),  (We  start  with  the  Fourier  integral  theorem  having  complex 
exponential  signs  opposite  to  87(a),  p 601). 
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Formally,  equation  (8)  does  not  converge  on  the  large  semicircle 
because  of  a factor  in  the  integremd  modulus  arising  from  (^*o  v^* 


We,  therefore,  actually  evaluate 


d+i» 


(‘‘v)-^  / '“pK)  ’("<>)] 


d-ioo 


Term  by  term,  the  right-hand  sides  of  equations  (5)  to  (7)  have  no 
poles  (see  sect.  5) , for  o y 0,  other  than  the  pole  introduced  into 


equations  (6)  and  (7)  by  o - iioe 


In  order  to  evaluate  equation  (9) , 


we  resort  to  the  ansatz  of  calculating  for  the  H,  evaluation,  am 

•p 


integrand  multiplied  by  l/(s  + e) , with  c infinitesimal,  then  subsequent 
differentiation.  The  operamd  is  a step  function  which  yields  a Dirac 
delta  function  (see  sect.  5) . 


s s 

If  we  use  the  pole  at  z = -a/e  for  the  evaluation  of  E and  E*, 

o o r 0 


the  factor  (a/v) (o/e^)  appears  in  our  expressions. 


We  can  evaluate  E®  and  E?  by  a similar  ansatz  used  for  H.  by 

IT  “ w 


extracting  f(s)/[e  (s  + a/e  )]  rather  than  f(s)  from  the  Laplace 
o o I*  \ 

inversion  integrand.  With  f(0)  = e'  lO)  = 0,  the  subsequent 


convolutions  yield  the  same  result  as  the  a/e  pole  evaluations  for 

o 


a a 

— - « 1 . 


^o 


(10) 


The  approximate  field  solutions  should  be  subject  to  the  usual 
skin-depth  restriction  we  give  below.  In  equation  (9) , it  is 
interesting  that  (o/e  ) (a/v)  enters  the  n=l  electric  field  via  the 


factor 


^1  + [(o/Eq)  (a/v) 


]Vio) 


in  which  the  difference  from  unity  is 


small,  even  for  (a/e  ) (a/v)  1. 

o 


Lerch's  Theorem" 


assures  the 


uniqueness  of  equation  (15)  below  in  the  electric  field  evaluations. 


*F.  B.  Hildebrand,  Advanced  Calculus  for  Engineers , Prentice  Hall 


(1949) . (We  start  with  the  Fourier  integral  theorem  having  complex 
exponential  signs  opposite  to  87  (a) , p 601) . 


J 


including  tfte  electron 


The  resulting  solutions  for  a step  in  f(ti)). 


charge  sign,  are 

H®’^®P(r,0,t)  « a D(t)  sin  0 

(11) 

^ A 

®r  **(r,0,t)  » cos  0 f expj"-  (t  - T)j  d(t)  dx 

(12) 

Eg  ^(r,0,t)  = ^ (f)^  sin  e f exp  (t  - TjD(T) 

dx 

(13) 

where 

dx" 

(14) 

• ^ .‘J’- (-.■)■ 

In  equation  (14) , is  an  integer  order  modified  Bessel  function.  The 
desired-  responses  are,  therefore, 
t 

“ J^(t  - t)  [^c^^  ■ dT,  j^(0)=  0 (15) 


where  g ^ represents  any  of  the 


functions  in  equations  (11)  to  (13). 


At  r=a,  the  H field  and  surface  current  from  equations  (11) , (14) , 
2md  (15)  is 


(t)  =!  sin  0 

Ko  (t)  = -aj  (t)  sin  0 
0 c 


. (15)' 
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3.  ANTENNA  COUPLING  ESTIMATES 


We  desire  to  estimate  the  short-circuit  base-current  contribution 
for  2U1  antenna  mounted  on  the  sphere  as  shown  in  figure  2,  due  to  the 
field  E®  . The  direct  coupling  of  E^  to  the  antenna  is  not  considered 
here. 


Figure  2.  Antenna  coupling  geometry. 


We  assume  in  these  calculations  that  the  skin  depth  in  the  air  is 
significantly  greater  than  twice  the  sphere  radius.  Thus,  skin  depth  is 

- <^>'’  ’<10®  t (seconds)  (meters)  .>  2a  (16) 

where  t is  the  time  of  interest. 


Ilk 


The  x"  integral  in  equation  (19>  is 


£ (0.2)  I j '/ [ c )exp  - 2e^  ~ ^ (u  -yjx^  - l) 


1(0.2)  j (y/u^  - l) 


/u2^ 


du  = 0.2  In 


. (4^)] 


where  the  0.2  comes  from  a constant  approximation  (an  overestimate)  to 
e Ij (x)  shown  in  Abramowitz  and  Stegun.^ 


Some  values  of  interest  are  shown  in  tcUsle  I. 


TABLE  I.  PARAMETERS  FOR  THE  BESSEL  INTEGRAL 


Thus,  in  the  total  integral,  the  factor  multiplying  the  exponential 
in  the  x"  integration  is 


[- 1 (^)]  * >/([4^])  " (4^) 

it  (H^)) 

2 

M.  Abramowitz  and  J.  A.  Stegun,  eds..  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards,  AMS  55  (1965). 
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He  may  use  equation  (21)  to  qualitatively  judge  if  the  second  tern  from 
equation  (19)  is  not  dominant  in  our  estimate. 

The  Laplace  transform  of  equation  (18)  is 


E®(s) 


% -[ir  (H^)]  [-(^>]  • 


(22) 


We  shall  treat  the  amtenna  as  a monopole  over  am  Infinite,  perfectly 
conducting  ground  plane.  To  model  the  antenna,  we  use  a transsiission 
line  analog  due  to  Stark,**  shown  in  figure  3.  The  tramsmission  line  is 
lossless  (o  = 0)  amd  the  chauracteristic  iiqpedamce  is 


(23) 


vAiere  r^  is  the  antenna  radius. 


**W’erner  J.  Stark,  Transient  Response  of  a hog-Reriodic  lintenna  Based 
on  Broad-band  Continuous-Wave  Measurements,  Harry  Diamond  Laboratories 
TR-1792  (April  1977). 
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The  Incident  field  for  the  transmission  line  from  equation  (22) 


_inc,  . 2 

E (s)  « ~ cos 

X e 

o o 

where  we  have  approximated 
typically  in  table'll. 


(24) 


o 

e]q>onentlal  as  shown 


TABLE  II.  EXPONENTIAL  APPROXIMATION  TO  (a  - 1 m) 


r 

(?)’ 

exp  1 

: I <-=>] 

1-2.08 

I-I.5 

!■' 

a 

1.0 

1.0 

1.0 

1 .0 

1.25a 

0.51 

0.59 

0.69 

0.78 

1. 5a 

0.30 

0.35 

0.^7 

0.61 

2.0a 

0.13 

0.12 

0.22 

0.37 

o a 

In  comparison,  -r——  ■ 0,02  . 

o 10  ^ 

The  transmission  line  response  (see  sect.  6)  to 

is  shown  in  figure  4. 
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I 

i 


a 


®xp^-  — (t  - ti)J  - exp[-  ^(t  - ti)J 

c o 
2L  " e 

o 


exp  - ^ t - exp  - ^(t  - ti)  exp^-  Ij;-  ti^ 

c 0 ‘ 

2L  ” e 

o 

The  constant  L is  given  in  equation  (31)  below, 
current  estimate  is 


ti  ^ t <.  2ti  . (26) 


The  short-circuit 


^sc^^’  “ I 3(t  rV)'  [^c^^  " 


(27) 


where  the  voltage  euid  current  polarities  (see  fig.  3)  eure 


3V 

3x 


-L 

o 


3i, 

at 


(28) 


where  is  line  inductance  per  vinit  length.  Equation  (27)  has  a 

Laplace  transform  identical  to  the  more  familiar  form  in  which  j -*■  E®, 

c 

-*•  I^.  Uniqueness  is  again  assured  by  Lerch's  Theorem. 


We  choose  the  following  test  pulse,  shown  in  figure  5 


Figure  5.  Incident  electric  field  E ^md  Compton  current  pulse 


From  equations  (26)  and  (27)  (see  sect.  6,  eq  59) 


* M 

W 1 HI 

30 

IS 

Equations  (30)  and  (32)  are  only  valid  for  the  following  restriction 
since  only  was  so  simplified. 


exp(-|)«  1 . 


(33) 


4.  SAMPLE  CALCULATION 

If  we  use  the  model  of  section  3 for  a sphere  near  a ground  plane, 
we  should  attempt  to  take  some  account  of  E^  in  the  presence  of  the 
ground.  Following  a suggestion  of  Bombardt,  we  force  by  adjusting 
^m'  ^r'  ^ representative  as  possible  of  some  prescribed 
near-ground  horizontal  electric  field. 


Now, 


e E^  + oE^ 
o z z 


-J 

c 


emd  with  our  test  pulse. 


exp 

U 

( 

'r 

1 

(34) 


(35) 


20 


(36] 


In  order  to  have  a peak,  we  must  fix  t^  and  select  ^ 

For  the  exaople,  we  choose 

z * 372  ohms, 
c 

a > 1 m, 

- =2.08  (table  II), 

A 


4 m,  emd 
20  ns. 


in  «diich  case  the  exponential  more  closely  approximates  ( — ^ . The 
i / \ 

electric  field  at  r=a  for  one  value  of  is  shown  in  figure  5. 
Hie  short-circuit  current  estimate  for  the  case  in  figure  5 is  shown  in 
figure  6. 


The  peak  value  of  I is 
sc 

|l  I (first  peak)  -v-  (1.9)1  (1)  |cos  0| 


(37) 


L«t  US  now  vary  the  exponential  attenuation,  - , while  keeping  the 

A 

remaining  parameter  set  fixed.  Results  are  shown  in  figures  7 and  8. 
The  peak  currents  vary  monotonically  upward  as  the  attenuation  is 
decreased. 

The  ideally  gamma-thin  sphere  model^  yields  a net  time-dependent 
charge  which  is  equal  to  the  tO'tal  charge  that  would  be  present  in  the 
volume  of  conducting  medium  excluded  by  the  sphere.  Our  approximate 
driving  electric  field  E*  yields  aero  net  charge  on  the  sphere. 

Botabardtf  Jr./  Electromagnetic  Rm»pon»m  of  a Sphtto  Ovot  a 
Groxmd  Plane  In  thm  Presence  of  Source  and  Conduction  Currents  in  the 
Air,  Bmrzy  Dimmond  Lnbotntoxioa  TR~1716  (Octobor  1975). 
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The  neglected  term  in  our  electric  field,  equation  (17) , is  not 
quite  negligible  for  higher  conductivity  values.  However,  at 
interesting  lower  conductivities  the  second  term  in  equation  (21)  has  a 
typical  value  smaller  than 


and  the  contribution 
first  approximation. 


(0.2)5  »>  , 

2 X 10“^^  3 X 10® 

in  equation  (19)  is  small 


0.15  ; (39) 

enough  to  neglect  in  a 


5.  CONTOUR  INTEGRATION 

Since  our  method  of  evaluating  the  inverse  Laplace  transforms  is 
approximate  and  somewhat  unconventional,  we  include  some  further  details 
in  this  section. 

The  Bromwich  contour  for  evaluation  of  equation  (9)  is  shown  in 
figure  9.  We  first  examine  the  Hankel  function  denominators  in 
equations  (5)  to  (7)  (term  by  term).  From  a Bessel  relation,  the 
varnishing  of  denominators  requires 

h,!l!  (Jta)  = n h^^dca),  h^^  (ka)  = j (ka)  + iy  (ka)  . (40) 

I*  j-  n n n n 

From  equation  (40)  we  obtain 

,(ka)j  (ka)  = 0 (41) 

n“  ± n n- 1 n 

in  contradiction  to  the  known  relation 


24 


the  Hankel  func 


tion  factors  in  our  residue  evaluations 


varnish.  For  convenience  we  work  with  the 


Hie  argument  of  the  Hauikel  functions  ccxnes 
from 


Bromwich  contour 


First  we  state  that  (term  by  term)  the  modulus  of  exp  (a/v)z  g(2) 
vanishes  on  CB  and  HA  where  g(z)  are  the  integrands  resulting  from  the 
fields  in  equations  (5)  to  (7) . The  integrand  moduli  also  vanish  on  AB 


Likewise,  the  integrand  moduli  also  vanish  on  CD  and  GH  as  r 
We  find  that  we  may  remove  as  a convolution  factor: 


the  known  transform  of  which®  yields  equation  (14) , and  still  obtain 
veuiishing  moduli  on  the  contour  extremes  for  the  remaining  Integremd 


factors . 


The  loop  integrals  OKFG  are  of  the  form 


-(v*)"  [.  . ] /^Legenare 

n factors 


1.3..(2n  + 1)  ^ (2n  + 3)  • 


[Fjj(|x|e"’^^)  - Fjj(|x|e^^)] 


where  F^^  is  a Hankel  function  factor — the  exponential  in  the  Hankel 

products  is  in  the  convolution  product  of  equation  (48).  We  find  that 

the  small  circle  integral  vanishes  as  p ^ 0 because  of  the  p factor  from 
16 

dz  in  z = pe  on  the  circle.  We  find  explicitly  that  the  difference  of 
the  functions  F in  equation  (50)  is  always  zero  (even  for  I k|  = 0) j 
hence,  the  loop  integrals  contribute  nothing  within  our  approximate 
problem  formulation. 


The  residue  evaluation  for  H including  the  ansatz  is 


Pea(H  (z)  — ^ 
V (|)  z + e / 


(48)  removed 


sin  e expj^-e(t  + ^ (50) 


For  G emd  d infinitesimal,  we  differentiate  equation  (50)  to  obtain 

= a sin  06  (t)  . 

Equation  (51)  convolved  with  equation  (14)  yields  equation  (11). 


N.  w.  McLachlan,  Complex  Variable  Theory  and  Transform  Calculus, 
Cambridge  (1963) . 
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Using  the  pole  at  the  residue  evaluations  for  the  electric 

fields  yield  ® 

■ ^(r)’ ’ft) 

. ^ (f)  009  9 exp^-  t)  = (E^)  (53) 

Convolving  equations  (52)  and  (53)  with  (14)  yields  (13)  and  (12), 
respectively. 

Utilizing  the  6 function  ansatz,  the  electric  field  evaluation 

yields  equation  (15)  with  j (35)  and  sin0D(T), 

/sl\^  ^ / 

2(— ) cos  0D(t)  for  Eg,  E^,  respectively. 


6.  TRANSMISSION  LINE 


The  transmission  line  solution  in  sections  3 and  4 is 
straightforward.  We  include  the  following  development  for  conpleteness. 


The  Laplace-domain  short-circuit  current  for  the  (lossless) 
ismission  line  of  fic 
toward  the  open  end)  is 


transmission  line  of  figure  3 (x^  -►  is  the  space  variable  increasing 


1 (s)  » — 

ro  z 


/*■<, 


o o 


, ^ / 28h\ 

c 1 + exp(- 

-r)  - "’p(-  f - «„])] 


where,  for  a step 


(54) 


From  a series  development  of  equation  (54)  ve  find 


I„(t)  = Ii(t)  + l2(t)  + IjCt) 


vrtiere 


i,,t,  o<  t<a 

c 

- (Ir)  - f)]  i . 

and  we  are  interested  only  in  the  case 

exp(-  ~)  « 1 . 

Now, 

exp(-^)x.f.ii)  4 -a|) 

and 

X,,t,.-^[..exp(-^^)][f^.^ 

Assuming  equation  (58),  the  response  (56)  is  shown  in  figure  4. 


7.  CONCLUDING  RENAKKS 

The  radial  scattered  electric  field  solution  obtained  in  section  3 
should  be  used  strictly  in  an  essentially  quasistatic  situation.  It 
appears,  however,  that  only  an  exact  analytic  solution  for  the 
Compton-driven  sphere  in  a dissipative  medium  will  provide  a satisfying 
resolution  of  this  problem. 
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